Observing different quantum trajectories in cavity QED 
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The experimental observation of quantum jumps is an example of single open quantum systems 
that, when monitored, evolve in terms of stochastic trajectories conditioned on measurements results. 
Here we present a proposal that allows the experimental observation of a much larger class of 
quantum trajectories in cavity QED systems. In particular, our scheme allows for the monitoring of 
engineered thermal baths that are crucial for recent proposals for probing entanglement decay and 
also for entanglement protection. The scheme relies on the interaction of a three-level atom and a 
cavity mode that interchangeably play the roles of system and probe. If the atom is detected the 
evolution of the cavity fields follows quantum trajectories and vice- versa. 



I. INTRODUCTION 

One of the most remarking properties of quantum me- 
chanics is that the state of a quantum system changes 
not only via the deterministic evolution given by the 
Schrodinger equation but also when it is measured. Al- 
though the system can be measured directly, in many 
situations it is probed by an auxiliary system (ancilla) 
that is then detected, providing information about the 
system of interest. This is a typical situation in cavity 
QED where atoms and light interact and the detection 
of either of them alters the state of the other. In the mi- 
crowave regime, for example, it is usually the atoms that 
behave as a probe for the field. This configuration has 
led to a number of experiments on fundamental aspects 
of quantum mechanics including the measurement of the 
decoherence of a cat state QND measurement of sin- 
gle photons [2], and the observation of quantum jumps 
of light [3]. 

This property becomes particularly interesting for 
quantum open systems where the environment that the 
system is coupled to plays the role of a natural bona fide 
ancilla. The time evolution of a single quantum system 
can therefore be probed by directly monitoring this reser- 
voir. This dynamics is going to be stochastic, as the evo- 
lution is conditioned on the measurements results, and 
can be mathematically described in terms of quantum 
trajectories, which are related to different physical ways 
to monitor the environment and to extract information 
about the system. 

The monitoring of the field leaving a damped cavity 
mode provides a good example of this stochastic dynam- 
ics. For instance, if a photodetector is used to collect the 
output of the cavity, the dynamics is better described in 
terms of quantum jumps where each click in the detector 
corresponds to lowering the number of photons inside the 
cavity by one. However, a completely different dynamics 
is found if the same propagating field is combined with a 



local oscillator in a beam splitter and a homodyne mea- 
surement is performed, in which case the time evolution 
of the damped cavity field is appropriately described by 
a continuous stochastic trajectory (quantum state dif- 
fusion). All these trajectories present very interesting 
scenarios that allow for the production of non-classical 
states of the cavity field as well as the protection of the 
purity of the cavity mode or of the entanglement shared 
by two or more modes. Note that if the lost photons are 
ignored, or the measurements averaged out, one obtains 
the usual master equation dynamics for the decoherence 
process that is responsible, for example, for rapidly turn- 
ing superpositions of coherent states into mixtures [1]. 
For this reason, each monitoring scheme described above 
is said to represent an unravelling of the master equation 
in terms of stochastic trajectories. 

Note that the unravellings above represent only a lim- 
ited set of the possibilities to measure the environment. 
In fact the master equation can be mathematically unrav- 
elled in infinite ways in terms of stochastic trajectories [5- 
6 and therefore one could envisage more general ways of 
monitoring the environment rather than only simple pho- 
todetection and homodyning. This freedom in defining 
unravellings has been recently explored in the context of 
entanglement decay and protection [7 15 , where mon- 
itoring schemes that combine different decay channels 
play a crucial role on the recovery of the mixed state en- 
tanglement dynamics in terms of trajectories jSj [9l fTT} [T2] 
and on the protection of entanglement conditioned on 
measurement outcomes [15j . It would be interesting then 
to propose realistic experimental scenarios where this va- 
riety of unravellings could be explored. 

While the monitoring of the field emitted by a leaky 
cavity mode or by a decaying atom is within current ex- 
perimental feasibility, this is not necessarily the most ef- 
ficient way and certainly not the most complete one to 
generate different unravellings. This scheme proves to 
be very limited in two very interesting situations: non- 
zero reservoir temperatures and unravellings that com- 
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bine detections and no-detections (jumps and no-jumps). 
The first case presents two problems: first how to deter- 
mine when a very large reservoir loses a single photon to 
the system and second how to distinguish a photon that 
comes from the system from one that already exists in 
the reservoir. In the second case, one faces the problem 
of physically superposing a click with a no-click in the de- 
tector. These limitations severely hinder the exploration 
of quantum trajectories in these systems. 

Even though the direct monitoring of a natural thermal 
environment remains a challenge, in cavity QED one can 
artificially engineer this and other reservoirs that could 
produce different kinds of quantum trajectories when 
measured. In the microwave regime, for example, beams 
of atoms crossing the cavity can be used to mimic a ther- 
mal dissipative reservoir for the cavity field [16-18 . The 
posterior detection of these atoms produces quantum tra- 
jectories for the cavity field, as analyzed in [17 both in 
terms of jumps and continuous diffusion processes. How- 
ever, these schemes also present limitations to the pro- 
duction of different unravellings. For example, while the 
combination of different channels can be easily accom- 
plished within the optical detection scenario by having 
the photons for each channel arriving at different ports 
of a beam splitter [9l |T5] , the situation for atomic detec- 
tion seems far more complex. In all previous proposals 
to engineer thermal reservoirs for cavities using atomic 
beams, decay and excitation channels correspond to two- 
level atoms entering the cavity either in the ground or 
in the excited state. A combined detection, which is 
utterly important for the applications proposed in [8]- 
(TTJ [15] , would then require some kind of interaction be- 
tween the atoms after they cross the cavity, which albeit 
not impossible, seems rather challenging. 

In this paper we show that simple modifications to 
these previous proposals lead to an alternative scheme to 
engineer a thermal reservoir in the context of microwave 
cavity QED that accommodates more general detections, 
hence enabling the simulation of a wide set of different 
classes of trajectories for a monitored atomic reservoir. 
In particular, we address the limitations raised in the two 
previous paragraphs. We then propose a complementary 
experimental setup where we invert the roles, letting the 
cavity field play the reservoir for the atoms. We rely on 
the Purcell effect to channel the atomic decays into the 
modes of a lossy cavity and then use the optical detec- 
tion proposed in [15] to monitor the reservoir in different 
ways. In this last case, we introduce an additional mea- 
surement possibility to implement the proposal in [15] 
that would circumvent the challenging task of collecting 
broadly emitted photons from atomic decay. 



II. ATOMS AS THE RESERVOIR 

We will start by briefly recalling the usual way to sim- 
ulate a thermal reservoir for a cavity field using a se- 
quence of two- level atoms as discussed in [TBI HZl IE]- 



The atoms, prepared in either the ground or excited 
states, cross a cavity where they interact resonantly via a 
Jaynes-Cummings Hamiltonian for a short enough time. 
By short time, we mean a time interval St much smaller 
than the inverse of the vacuum Rabi frequency of the sys- 
tem 1/A such that the probability of an atom making a 
transition is small. Under this model, one can then show 
that if the atoms are not detected after interacting with 
the cavity, the cavity field (of mode annihilation operator 
a) will evolve according to a master equation of the form 

/) = 7_P[a]p + 7+D[at]p, (1) 

where T>{c)p = cpc^ — l/2{c^cp + pc'^'c), and the relation- 
ship between the rates 7_ and 7+ is given by the ratio 
between the flux of atoms initially prepared in the ground 
ivg) and excited state (rg). A thermal bath is obtained 
by setting this ratio to be Ve/vg = n/(l + n), where n is 
the average number of photons in the modes with energy 
hjOc in the reservoir. 

If one now considers that the atoms are detected after 
the interaction, then the evolution of the cavity field is 
conditioned on the result of this measurement and bet- 
ter described by quantum trajectories. In this case, if an 
atom enters in the ground (excited) state and is detected 
in the excited (ground) state then the cavity field will be 
modified by a jump operator J_ = a (J+ = a"*") corre- 
sponding to the annihilation (creation) of a photon in the 
cavity. Note that these jumps will be rare since, under 
our assumption of small interaction time, the atoms will 
most probably remain in their original state, in which 
case the cavity field evolution will correspond to the ap- 
plication of a no-jump operator Jq = (1 — y If 
the state of the atom is known before and after the cav- 
ity then a quantum jump unravelling corresponding to 
the operators {Jq, J-, J+} will be produced for the evo- 
lution of the cavity field [17 . Note, as mentioned before, 
that in order to produce any other jump-like unravelling 
which necessarily combines the operators J_ and J+, the 
atoms would need to interact after going through the 
cavity which is a rather difficult experimental task. 

We now present an alternative scheme to engineer a 
thermal reservoir that naturally allows more general de- 
tections. The central idea is to encapsulate in a single 
atom the interactions with the cavity that correspond to 
the evolution with J_ and J+. For that we will use three- 
level atoms with levels organised in a cascade, as shown 
in Fig. [1] and selected in such a way that the energy 
difference between the intermediate level (|e)) and the 
lower and upper levels (respectively \g) and \i)) are close 
enough to the energy of the chosen cavity mode so that 
these transitions can be tuned into resonance through 
the application of external fields, but only one at a time. 
The atoms are prepared in the |e) state before enter- 
ing the cavity and in the first stage, the external field 
shifts the |e) level such that the transition \g) \e) be- 
comes resonant with the cavity frequency cJc- The sys- 
tem then evolves under the Jaynes-Cummings Hamilto- 
nian Hi = hXi{\e){g\a + \g){e\a^) for a short time Sti. 
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Up to second order on X^Sti the initial atom- field state 
\^{t)) = |e)|i>) will be transformed to 

+ 5h)) = (|e) - (Ai«i)!aat|e) - iX,5h a^g)) |$), 

(2) 

where Ai is the couphng constant, |$) the initial field 
state and the tilde indicates that the state is not nor- 
malised. If, now, the other transition |e) \i) is tuned 
to resonance by changing the external field, the Hamil- 
tonian will be H2 = hX2{\i){e\a + \e){i\a^) and, again, 
making the short time expansion will lead to 

\^{t + Sti + St2)) = (|e) - iXiSti a^g) - 1X2^12 a\i) 

^ ^ a^a\e) \ \^). (3) 

Considering now a time interval At that is large enough 
so that n atoms (n = rAt ^ 1) cross the cavity and yet 
the probability of occurrence of a quantum jump remains 
very small, one can recover the master equation result 
by tracing out the atomic system. This becomes more 
evident if one identifies the rates as 7+ = r{Xi5ti)'^ and 
7_ = r(A2(5t2)^, where r is the number of atoms entering 
the cavity per unit of time. Therefore, the same atom 
can emulate at the same time both the dissipative and 
the excitation reservoirs. 
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FIG. 1: Atoms with the level scheme shown on the left cross 
the cavity C before being detected. Inside the cavity, the 
application of external fields change the atomic levels in two 
stages: first the cavity mode becomes resonant with the e — g 
transition, then with e — i transition. Rotations between levels 
i and g after the interaction with the cavity and before atomic 
detection produces different unravellings of the dynamics in 
terms of quantum jumps. 

If instead of tracing out the atoms one measures their 
levels, then the evolution of the state of the cavity field 
obeys a quantum trajectory given by the no-jump opera- 
tor Jo = (1 — -^^^aa^ — ^~^^ a^a) if the atom is detected 
in level |e) or the jump J_ = y^7_ At a (J+ = y^j^~Ki a^) 
if the atom is detected in level \i) (|^)), that corresponds 
to the usual monitoring discussed previously. 

The advantage of encoding both jumps in the same 
atom is that it allows us to explore different unravel- 
lings by introducing simple modifications in the detection 
scheme. If one makes a unitary transformation between 



levels \i) and l^f) after the atoms cross the cavity, for ex- 
ample, then a detection of each state will correspond to 
a different unravelling in which the jump operators are 
combinations of a and . For instance, in the partic- 
ular case that 7_ = 7+, a 7r/2 pulse between levels \g) 
and \i) produces jump operators proportional to the field 
quadratures a + e'^^a^ and a — e~'^^a^ where the phase 
(j) is given by the chosen rotation. In fact, since one can 
apply unitary operations also involving level |e) then it is 
clear that any unravelling of the type U J, where [/ is a 
unitary matrix acting on the vectorial space of the jumps 
and J is the vector J = {Jq^ J-^^ J-)^ can be produced. 
Note that the unitary U is obtained by individually ro- 
tating the atoms after they have already interacted with 
the cavity mode therefore preserving the master equation 
structure when the atoms are ignored. 

In the scheme described above, the detection of the 
atom modifies the quadratures of the field affecting, in 
principle, all Fock states. For example, if the field is ini- 
tially constrained to Fock states {|0), |1)}, the detection 
of the atom in its middle (|e)) or highest levels will 
preserve the initial subspace, but the detection of the 
atom in its ground state will expand the effective sub- 
space of the cavity mode to include Fock state |2). This 
can be avoided by choosing a selective interaction [21] 
for the |e) \g) transition in which case the creation of 
an extra photon in the cavity can be set to occur only if 
the cavity is empty, hence preserving the initial effective 
subspace and thus fully mimicking a spin-half system for 
the cavity field [22j . 



III. CAVITY FIELDS AS RESERVOIRS 

Inverting the roles of the previous session, now, the 
three level atom is stationary and interacts with two or- 
thogonal cavity modes such that the \i) \e) (|e) \g)) 
transition generates photons circularly polarized to the 
Right (Left) associated to the annihilation operator 
{aL) of the cavity mode. We assume a very lossy cavity 
(of decay rate n) made of a nearly perfect mirror and a 
semi-transparent one (see Fig.|2|. If the Rabi frequencies 
of each transition are large enough then the Purcell effect 
will channel the atomic decay in the cavity modes, and 
the cavity asymmetry will guarantee that the photons 
will always leak in a well defined direction. 

Each atomic transition has to be coupled to its re- 
spective cavity mode with different coupling constants 
Xie and Ae^, and a 7r-polarized classical field of strength 
Vt has to drive the \g) \i) transition as shown in 
Fig. [2] The relation between the involved rates should 
he Xie > ^eg > ^ ^ 7, T, where 7 and F 

are the natural atomic linewidths. For the scheme to 
work, much greater typically reads (z^, A^g, Ag^, 7, F) = 
(1000,100,5,1,0.01,0.01) where this particular choice 
has already been made to fit the analysis that follows. 

The overall master equation that describes the dynam- 
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FIG. 2: An atom with the level scheme shown on the left is 
placed inside a cavity and is coupled to two orthogonal modes 
of a damped cavity (right panel). The atom-cavity couplings 
{Xge and Xie) and the cavity decay (k,) will induce an effective 
atomic decay via the Purcell effect. Photons from the right 
and left polarized modes leaving the cavity will then corre- 
spond to |e) \g) and \i) \e) transitions, respectively. 
The corresponding effective decay rates jeg and jie are as- 
sumed to be much larger than the natural spontaneous decay 
rates 7 and F. 

ics of the whole system is given by: 

p=-'-[H,p]+^V[a^y^]p + rV[a^% + KV[a], (4) 
where 

— = uje\e){e\ ^uji\i){i\ ^URO^^aR ^ULa^aL 

-^i{cj^;'^e'^^'-a^_l'^e-'^^'), (5) 

with cr!f and cr!f? being the lowering and raising opera- 
tors for the ij transition. Being hz Xge there is no 
Rabi oscillation between the corresponding atomic tran- 
sition and the cavity mode which acts just as a channel 
for the atomic decay. In other words, as soon as a pho- 
ton is transferred from the atom to the cavity mode it 
leaks through the semi-transparent mirror. The effective 

atomic decay rate is then given by 7_ = jeg = as 

long as 7 <C Much in the same way, the effec- 

tive decay rate for the \i) \e) transition is given by 

Tie = Now, if the external driving field Q is still 

much smaller than this effective decay rate, then level \i) 
can be adiabatically eliminated and the combined effect 
of the classical field and the strong decay rate will be 
to generate an effective incoherent pump between levels 
1^) and |e) of rate 7+ = ^ = In this way, the 

scheme produces a decay and an excitation reservoir for 
the qubit composed of levels e and g with rates 7_ and 
7+, respectively. In particular, when the effective rates 

7± are the same, i.e. = one obtains an infinite 

ie 

temperature environment. 

Note that there is still another degree of freedom avail- 
able if the conditions are not perfectly matched that is 
the eventual detuning between cavity modes and atomic 
transitions. For example, let us say that Xge ~ Xei- If 



Age = cje — > 0, then the effective decay rate is cor- 

rected by 7_ = ,(i+4^|^/,2) 

Under the above approximations and assuming that 
the effective rates are much larger than the natural 
atomic linewidths (7^ 7, T), the unconditioned dy- 
namics of the (two-level) atom can be finally written as: 

p^7_P[a_]p + 7+D[a+]p, (6) 

where the superscript {eg) has been dropped for simplic- 
ity. If the photons leaving the cavity are detected, then 
the atomic dynamics will follow a stochastic dynamics 
conditioned on the measurement outcomes as described 
by the quantum trajectory approach. Equivalent ly to the 
previous session, different detection techniques will lead 
to different unravellings. For example, if a A/4 plate is 
placed right after the leaking mirror, then each circular 
polarization is converted into orthogonal linear polariza- 
tions, aR ^ an {cll ^ ay). These photons are then sent 
into a polarized beam splitter that separates the linear 
polarizations into two different propagation modes where 
detectors are placed to collect them. In this way, the 
click in the an {ay) path identifies a transition \g) \e) 
(|e) \g)) which corresponds to the "usual" photodetec- 
tion unravelling with the jumps J+ = y^^7+^<^+^^ and 
J_ = Y^^TSai^^^ applied to the atomic state. If no 
photon is detected in the dt interval, then the no-jump 
Jo = (1 — -^^^(Jee — ~^^^^gg) Operator is applied to the 
system. 

A more interesting situation occurs if the A/4 plate is 
taken out of the setup. In this case, the photons propa- 
gating after the PBS will correspond to linear combina- 
tions of the photons leaking from the atom, implement- 
ing an unravelling with the jumps ax = ai^^^ + a^^^ and 

ay = cr^^^ — cr^^'^ (similar to the quadrature unravelling 
described in the previous section) . We have shown in |T5] 
that if two qubit s initially share an entangled state, then 
this kind of monitoring performed on both subsystems 
can preserve the entanglement in the system. The above 
setup is therefore a way to implement the entanglement 
protection ideas proposed in [15 . 

Before concluding, let us just remark that the scheme 
here described for cavity QED could also be applied to 
other experiments involving three level systems and har- 
monic oscillators. Natural candidates where quantum 
jumps have been recently observed or proposed are su- 
perconducting qubits [23l|24] and nanoresonators [251126]. 

Finally, in conclusion, we have proposed two experi- 
mental ways to produce a wide range of unravellings of a 
master equation evolution in cavity QED systems. In the 
first case, three- level atoms are used to simulate a ther- 
mal dissipative reservoir and the preparation and poste- 
rior detection of these atoms in different basis produce 
unravellings that correspond to different combinations of 
the natural jump operators and even between these and 
the no-jump one. This is a particular case in which "no- 
click" and "click" can be physically combined to generate 
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a new class of unravellings. Later, we also show how to 
invert the roles and use the Purcell effect to channel the 
detection of the spontaneous decay and of the incoherent 
pump (or even their combination) for driven three-level 
atoms. This work adds to the existing class of proposed 
cavity QED experiments, expanding the range of quan- 
tum effects that can be explored in this experimentally 
successful system. 
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